00 

o 
o 

(N 



Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 10 July 2008 (MN LM^ style file v2.2) 

Deciphering the Pulsations of G 29-38 with Optical 
Time Series Spectroscopy 

Susan E. Thompson^'^*, M. H. van Kerkwijk^, and J. C. Clemens^ 

^ Dept. of Physics and Astronomy, University of Delaware, 217 Sharp Lab, Newark, DE, 19716 USA 
'^Dept. of Physics, Colorado College, Colorado Springs, CO 80903 USA 

^Dept. of Astronomy and Astrophysics, University of Toronto, 50 Saint George St., Toronto, ON M5S 3H4, Canada 
Dept. of Physics and Astronomy, University of North Carolina- Chapel Hill, Chapel Hill, NC 27599 USA 

6 

^ ' ABSTRACT 

We present optical time series spectroscopy of the pulsating white dwarf star 
QQ ' G 29-38 taken at the Very Large Telescope (VLT). By measuring the varia- 

vQ I tions in brightness, Doppler shift, and line shape of each spectrum, we explore 

the physics of pulsation and measure the spherical degree (i) of each stellar 

c^ ■ 

f~^ ■ pulsation mode. We measure the physical motion of the g-modes correlated 

oo ■ 

o 



with the brightness variations for three of the eight pulsation modes in this 
data set. The varying line shape reveals the spherical degree of the pulsations, 

•'r-t . 

^^ , an important quantity for properly modeling the interior of the star with as- 

j^ ■ teroseismology. Performing fits to the H/3, H7, and H5 lines, we quantify the 

changing shape of the line and compare them to models and previous time 
series spectroscopy of G 29-38. These VLT data confirm several £ identifica- 
tions and add four new values, including an additional £—2 and a possible 
£=4. In total from both sets of spectroscopy of G 29-38, eleven modes now 
have known spherical degrees. 

Key words: stars: oscillations - stars: white dwarfs - stars: individual : 
G 29-38 
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1 INTRODUCTION 

DAVs, or ZZ Cetis, are variable white dwarf stars with hydrogen atmospheres. They are the 
coolest known white dwarf pulsators and reside in an instability strip near 12,000 K. Their 
pulsations are multi-periodic with amplitudes as high as 5 per cent and periods between 
70 s and 1100 s. In the past few ye ars the number of known DAVs has increased from a few 



dozen to more than a hundred (see 



2005 



Mukadam et al 



Voss et al. 



12007 



Castanheira et al. 



2006 



MuUallv et al. 



20041 ). The pulsators at the hot end of the instability strip show lower 



amplitudes and more stable pulsations, while the others show more complex pulsations 
including variability of the excited modes, and an abundance of harmonic and combination 
modes. 



The high gravity (log(g)~8 cgs) of the DAVs favor gravity- mode type pulsations (JRobinson et al 



19821 ). where motion along the surface of the star produces regions with varying temper- 
atures. To describe the spatial distribution of these regions, each of the pulsations on the 
star are characterized by a spherical harmonic (Y^ „i) and the radial order, n. The spherical 
degree, i, and azimuthal order, m, describe the number of nodal lines across the surface of 
the star, and directly affect observations of the pulsations. 

The pulsation's radial component allows asteroseismological models to determine the in- 
ternal structure of the star. The DAVs have evaded unambiguous asteroseismological mod- 
eling in part because of the difficulty in measuring the spherical degree of their pulsation 
modes. Usually a value of i=l is assumed in these models. However, it is possible to fit 



modes with h igher va 



For e xample, 



Bradley 



u es oi £ and this can significantly change the parame ters of the model 



20041 ). instead of i=2 (jCastanheira et al. 



(120061) showed that usin g the identification of £=4 (JThompson et al. 



20041 ) for one mode, changes the modeled depth of 



the Hydrogen layer in G 185-32 by a factor of 100. 

One method to determine i independent of asteroseismological modeling was introduced 



by 



Robinson et al. 



(Il995l ). They showed that i can be obtained from time series optical and 
UV spectroscopy. The effects of limb darkening increase in the blue wavelengths and decrease 
through the spectral lines, changing the fraction of integrated light coming from the limb. 
As such, t he pulsation amplitude varies with bo th wavelength and the spherical degree of 
the mode. IClemens. van Kerkwijk. fc Wul (120001 ) (hereafter COO) measured chromatic am- 
plitudes (pulsation amplitude at each wavelength bin) across the Hydrogen Balmer lines 
and demonstrated that spectral variations can indeed be used to identify spherical degree. 
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^he same method was applied to the DAVs HS 0507+0434B, G 117-B 15A, and G 185-3 2 



(JKotak et al. 



2002 



2004 



Thompson et aLll2004h and the DBV GD 358 (jKotak et all 120031 ). 



Generally, the method worked well, although for one r node - the 272 s rn ode in G 117-B15A 
- the chromatic amplitudes are extremely puzzling (jKotak et al.ll2004l ). Also, they found 
that the chromatic amplitudes provided unambiguous identifications only for the strongest 
modes; for w eaker ones, noise ob scured the expected visible differences. In order to improve 



the precision. 



Kotak et al 



(120031 ) attempted to use the slope and curvature in selected wave- 
length ranges of the chromatic amplitudes as measures of s pherical degree, Av i th som ewhat 
limited success. A different approach was demonstrated by [Thompson et al.l (120041 ) where 
they analysed fits to the spectral lines rather than the observed fiuxes in each spectral bin. 
The chromatic amplitudes created from these fits did allow for mode identifications from 
lower signal-to-noise data. 

G 29-38 is one of the brightest DAVs with complex pulsations and some of the largest 
amplitudes, making it a prime candidate for time series spectroscopy. It has an abundance of 
modes, but they are not all excited simultan eously and many observation runs are necessary 



to measure them all. 



Kleinman et al 



(119981 ) compiled 10 years of photometric observations 
and found a possible series of 1=1 modes on G 29-38. Time series spectroscopy taken at the 
Keck Observatory with the LRIS spectrograph provided chromatic amplitudes of G 29-38 



that d irectly detected an £=2 mode at 776 s ( COO) and gave strong evidence 



920 s (IKotak et al. 



2OO2I ). Using the same data. Ivan Kerkwijk. Clemens, fc Wul (I2OOOI ) (here 



'or another at 



after VKOO) measured the periodic line-of-sight velocities associated with the longitudinal 
motion of the gas. For the 776 s mode they found relatively large motion as compared to 
the measured fiux amplitude, again indicating that the mode is £=2 ( Dziembowskilll982l). 



The measurements of the pulsation velocities were confirmed by [Thompson et al. 



(12003h 



using a series of high resolution spectra of G 29-38 and provided observational evidence 



regar ding the convective driving theory proposed for DAVs (jBrickhill 



1983 



Goldreich&Wu 



I999I ). The velocities also hint at the nature of the many combination and harmonic modes 



present on these pulsators. These combinations and harmonics most likely result from r ion- 
linea r mixing at the sur 



2001 



Montgomery 



ace of the star and are not real modes that probe the interior (IWu 



20051 ). VKOO reinforced this idea by finding no velocity signatur e asso - 
ciated with the combinations and harmonics on G 29-38. However, [Thompson et al.l (120031 ) 
reported one combination mode with a significant velocity detection. 

In this paper we show the results of time series spectroscopy of G 29-38 taken at the 
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VLT (Very Large Telescope) in 1999. We present the light and velocity curves measured 
from the series of spectra. We fit 13 modes in the light curve and find variations at the 
same period in the velocity curve for the largest amplitude pulsations. We then perform 
fits to the Balmer lines to quantify the changing line shape. Since we are using fits to the 
observed spectral lines, we can present the line shape variations in terms of the fitted values: 
the normalized equivalent width variations of the Gaussian and Lorentzian functions. We 
perform the same analysis on the 1996 Keck spectral series presented by VKOO and COO as a 
check of our method and recover their results of an £=2 mode among a series of i=l modes. 
Application of the same method to our VLT data adds four additional mode identifications 
to the G 29-38 pulsation spectrum. 



2 OBSERVATIONS 



We present a 6.14 hour spectral series 



604 spectral image s) of G29-38 from the VLT, 



observed with the FORSl spectrograph (JMoehler et al 



19951 ) on August 27, 1999 starting 



at 03:00:19.67 UT. The average exposure time on each image is 16.0 s with a readout time 
of 16.0 s. We perform a standard reduction on each image using routines in IRAF. We 
subtract the bias measured from the overscan, divide by a normalized average flat taken at 
the beginning of the night, and perform wavelength calibration from arc lamps taken at the 
end of the night. The average flat field was insufficient to remove all the effects of erroneous 
pixels, thus small dips and peaks appear in the average spectrum. Overall, these are smaller 
than the average noise in the individual spectrum. Also since we rely on relative fluxes, poor 
fiat-fielding does not significantly change our results. No spectro-photometric standard was 
available to flux calibrate the spectra. Each extracted spectrum spans 3300 A to 5700 A, has 
a dispersion of 2.35 A per binned pixel, and a signal-to-noise of 80. Because light entered 
through a wide slit, the resolution is set by the seeing, yielding ~8 A. 

Figure 1 depicts an example of an individual spectrum and the average of the 604 spectral 
images. Besides the Hydrogen lin es ranging from H/5 to H9, the Ca K line is evident at 3933 A. 



von Hippel fc Thompson! (120071 ) use all known calcium line equivalent width measurements. 



including those presented here and in VK OO, to show that the abundance of Ca on G 29-38 



varies by as much as 70 per cent, though 



Debes &: Lopez- Morales! ( 12008! ) has since seen no 



variations in the Ca of G 29-38. Measurements of changing Ca abundance will allow us to 
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Figure 1. Average and individual spectra taken with the FORSl spectrograph. All small dips in the average spectrum are a 
result of poor flat-fielding, except for the Ca II line at 3933 Aand the Mg I line at 4481 A. 



;he WD atmosphere and t 



constrain the diffusion rate of calcium in 

from the debris disk surrounding G 29-38 (IZuckerman &: Becklin 



le variable accretion 



1987 



Reach et al 



20051) 



3 PERIODICITIES 

Time series spectroscopy of pulsating white dwarfs provides the opportunity to measure both 
the variations in brightness of the star and the radial velocity of the gas along the surface 
of the star. We find periodic variations in our measurements of both the light and velocity 
curves. We create a light curve by averaging the counts collected in the continuum red-ward 
of the H/3 line (5000-5647 A). The average noise level of the light curve is 0.89 mma (0.089 
per cent), measured from the high frequency end of the Fourier transform (FT). The velocity 
curve is created from the weighted average of the Doppler shifts, found by individually fitting 
the H/?-He lines with a combination of a Gaussian and Lorentzian function fixed to have 
the same central wavelength. The average noise level of the velocity FT is 0.56 kms~^. See 
Figure 2. 

We fit and remove individual sine waves to the light curve starting with the largest ampli- 
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tude mode, the order is depicted in Table 1. We then fit modes that represent combinations 
of these larger modes. Overall, we find eight independent frequencies and five combination 
frequencies in the light curve. The final fit in Table 1 represents fitting all 13 modes si- 
multaneo usly. The errors represent analytical est imates of the uncertainties calculated by 



Period04 (JLenz &: Bregerll2005l : iBreger et al.l 119991 ). a program that performs FTs and sine- 



wave fits to a given light curve. The reduced chi-square value for the fit to the light curve is 
2.4, indicating that our 13 modes do not represent the entire signal in the light curve given 
the size of the error in each light curve measurement. The remaining signal is dominated by 
low-frequency signals that are typical of variable seeing and are probably not inherent to the 
star. The nonlinear light curve of G 29-38 also suggests that combination modes below the 
noise level are present proving further signal that we did not remove. The errors presented 
in Table 1 represent this reduced chi-squared value. 

The strongest periodicity of the velocity curve corresponds to the strongest mode (Fl) 
found in the light curve, giving us confidence that we have measured the motion of the gas 
on the surface. The light curve measures the period of each mode better than the velocity 
curve because of its lower noise level. We fit the velocity curve with sine waves, fixing the 
frequencies to the thirteen modes measured in the light curve. The reduced chi-square for 
the fit of the 13 frequencies to the velocity curve is 1.6. The remaining signal in the velocity 
FT is at low frequencies and probably result from stellar wander in the wide slit. The lower 
reduced chi-square of the velocity fit as compared to the light curve fit indicates that the 
velocity curve is more linear, i.e. no undetected combination frequencies. The amplitudes 
and frequencies found by fitting the velocity and light curves at the same frequencies are 
marked on the FTs in Figure 2. 



To determine the significance of a mode we use the false alarm probability (lKeplerlll993l ). 
In the light curve, to have less than 0.1 per cent chance of a peak being produced by noise 
alone it must have an amplitude of at least 3.2 mma. For the combination modes, where 
we know the frequency lies at the combination of the two dominant modes, the peak need 
only have an amplitude of 2.3 mma. Table 1 only includes peaks in the FT with amplitudes 
greater than these values. For the velocity FT, peaks need an amplitude of 2.04 kms~^ to be 
detected independent of the light curve with the same 0.1 per cent confidence level. However, 
since the light curve fixes the frequency, for a velocity mode to have less than 0.1 per cent 
chance of being due to noise it simply needs an amplitude greater than 1.47 kms~^. Only 
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Figure 2. The Fourier transforms of the hght (top) and velocity (bottom) curves of G 29-38. The modes found in the data 
are labelled at the top. Sine-wave fits to each curve at these frequencies are marked on the FT. Pluses represent independent 
modes and crosses represent combination frequencies. The residuals after the fit are shown as reflections beneath each FT. 



the modes Fl, F2 and the combination of F1+F2 have velocity amphtudes greater than this 
hmit. 

We compare the observed velocity and light amplitudes and phases with the values 
Rv =(A^,/ci;)/A/, relative amplitude, and A$ = $„ — $/■, relativ e phase in Table [TJ The 
measured R.^ values have been shown by [Thompson et al.l (120031) and VKOO to generally 



agree with the convective driving theories of 



Goldreich fc Wul (Il999l ). Our measurements of 



G 29-38 agree with the previous measurements with the Keck telescope. One mode present 
in all three data sets (615 s) co nsistently has the s ame r elative amplitude: R„ is 17±3, 12±4, 
and 18±2 Mm/rad for VKOO, 



Thompson et aL 



(120031) and this pape r respectively. These 



measurements are all near to the prediction by lWu fc GoldreichI (119991 ) of 16 Mm/rad for an 
800 s mode on a similar temperature star. The relative phases, A$, for those modes with a 
significant velocity dete ction all show the expec ted lag of the brightness behind the motion 



of the gas, predicted by 



Wu fc GoldreichI (119991 ) to be a 55° lag for an 800 s mode. 



With the observation of several combination modes, we compare the combination ampli- 
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Table 1. The modes measured in the Ught and velocity curves of G 29-38. The fits to the flux time series use the function AfCos{2TTft — 
^f) and the velocity time series fits use the function AyCos{2nft — ^v)- Comparison of the flux to the velocity pulsations are R^ = 
{Ay/Lij)/Af and A<1? = <l?i, — <I>j. The comparison of the combination modes to their parent modes are given as Re = A\ ■' /riijAKA''. 



and A-I-c = *)^+^ 



-*P 



P 



/ 
(/.Hz) 



(mma) 



n 



Av ^v 

(kms-1) {°) 



Rv 

/ Mrn\ 



A$ 

n 



Re 



A*c 

n 



Fl 

F2 

F3 

F4 

F5 

F6 

F7 

F8 

2F1 

2F2 

2F4 

F1+F4 

F1+F2 



615 
811 
835 
353 
937 
681 
283 
773 
308 
405 
177 
224 
350 



1625. 
1233. 
1197. 
2829. 
1067. 
1469. 
3534. 
1293. 
3250. 
2466. 
5658. 
4454. 
2858. 



0±0.9 

4 ±1.0 

2 ±1.0 
1±2.1 
6 ±2.1 

3 ±3.0 

5 ±5.8 
9 ±4.3 
1±0.0 
8 ±0.0 
2 ±0.0 
1±0.0 

4 ±0.0 



27.7 ±1.0 
22.9 ±1.0 

13.4 ±1.0 
10.7 ±1.0 

10.5 ±1.0 
8.3 ± 1.0 
4.1 ± 1.0 
5.9 ±1.0 
3.0 ±1.0 

7.5 ±1.0 
2.3 ±1.0 
3.7 ±1.0 

2.6 ±1.0 



226 ±2 
189 ±2 
192 ±4 
82 ±5 
290 ±5 
274 ±6 
347 ± 13 
238 ± 10 
76 ±19 
356 ±7 
110 ±23 
283 ± 15 
70 ±24 



5.2 ±0.5 
2.1 ±0.5 
0.5 ±0.5 
0.3 ±0.5 

1.3 ±0.5 

1.1 ±0.5 

1.2 ±0.5 
0.3 ±0.5 

0.48 ±0.5 
0.37 ±0.5 
0.72 ±0.5 
0.37 ±0.5 
1.74 ±0.5 



286 ±5 
225 ± 13 
102 ± 55 
52 ±54 
33 ±20 
341 ± 24 
37 ±20 
23 ±90 
180 ± 55 
83 ±70 
347 ± 36 
114 ±71 
87 ±15 



18.3 ±1.7 
11.8 ±2.8 
4.9 ±4.8 

1.6 ±2.5 
18.8 ±6.8 
14.3 ±6.2 
13.8 ±6.1 
6.3 ±9.6 

7.7 ±7.8 
3.2 ±4.0 
8.6 ±6.5 
3.6 ±4.5 

37.8 ±17.2 



60 ±5 
36 ±13 



103 ± 21 
67 ±25 
51 ±24 



17 ±28 



3.9 ±1.3 

14.1 ±2.1 
20.6 ±9.2 
6.2 ±1.8 
2.0 ±0.8 



-16 ±19 
-22 ±8 
-54 ± 24 
-25 ± 16 
15 ±24 



WJ (J200ll). Specifically, 



tudes and phases with their parent modes according to VKOO and 

we show the combination amplitude ratio, Re = A^j^-' /riijAlA-'f, and the combination phase 



$}-^- 



difference, A$n = ^Y'' — ($} + $}), in Table 1. The Re values for the first harmonics of 



our Fl and F2 agree within the errors with the same modes presented by VKOO on G 29-38. 
Since combination modes do not suffer from cancellation over the surface of the star in the 



same manner a s their parent modes, the value of Re can be an indication of £ (JYeates et al 



2005 



Wu 



20011 ) ■ However, these values are dependent on the azimuthal order, m, and the 



inclination of the star. Since we have not resolved rotational splitting of our modes, use of 
combination modes to determine i is very uncertain. We only note that the first harmonic 
of F4 shows a larger value than the other harmonics. 

One curious feature in these time series is the modestly significant peak in the velocity 
FT at 2858//Hz that matches the combination frequency FH-F2. While the combination of 
F1+F2 lies within 30/iHz of F4, it is impossible to satisfactorily fit the observed velocity 
peak with the same frequency as F4. While the 1.74km s~^amplitude of this velocity curve is 
above the significance level, this criterion only applies if we can expect to find combination 
frequencies in the velocity curve. Otherwise the false alarm probability is 1/26, i.e. a 3.8 per 
centchance of the velocity peak being due to noise. 

Most observations and theory insi st that combination modes should have no motion 



corresponding to the flux variations ( jWu 



2001 



van Kerkwijk et al.l l2000l ) . Combinations 



are believed to occur from the changing 



delaying and distorting the pulsations (JMontgomery 



depth of the oute r convectio n zone, attenuating 



2005 



Wu 



200 ll ). Hence we observe 



nonlinear pulse shapes in the light curve and combination modes in the brightness FT. 
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Since the combinations occur from nonlinear mixing of the temperature v ariations, they do 



not re quire gas motions to be seen as a peak in the hght FT. However, 



Thompson et aL 



(120031) previously found one instance of a significant velocity variation on G 29-38 that 
matches the sum of two large modes. So, despite the theoretical predictions, observations 
tell us that we can expect to find velocity combination and we conclude that F1+F2 is the 
second observed combination mode demonstrating an associated periodic velocity variation. 



4 MODE IDENTIFICATION OF G 29-38 

Due to differences in limb darkening across the spectral lines, the hydrogen spectral line 
shapes change according to the distribution of the pulsation across the surface of the star. 
As a result, chromatic amplitudes, wavelength dependent amplitude variations across the 
Balmer lines, have characteristic shapes that reveal L This technique cann ot reveal m be- 



cause of the purely spherical harmonic treatment of the brightness variations (jRobinson et al 



19951 ). Also, the motion of the gas, which is variable over the surface of the star and depen- 
dent on m, is too small to observe except for the average change of the line center as shown 
in the previous section. 

Using the technique of chromatic amplitudes, COO show that one mode has £=2 while all 
others match £=1 on G 29-38. We present the chromatic amplitudes from the VLT data of 
the six largest pulsation modes in Figure [31 As seen across the H/? through H5 lines, our data 
demonstrate the typical central peaks with side dips predicted by the atmospheric models 
for low values of L However, the noise, especially in the lower amplitude modes, prevent 
us from making specific mode identifications when we compare them to the models. One 
obvious feature is that the 353 s mode, F4, appears to have shallower dips on either side of 
its peaks, indicating that F4 is a different spherical degree than Fl or F2. F5 also seems to 
have somewhat shallower dips, though the increased noise makes the difference much less 



significant. 

Since the introduction of chromatic amplitudes, fitting spe ctral lines has proven 



to be 



useful in extracting spherical degree from lower quality data (iThompson et al.ll2004l ). As 
a result, we fit the spectral lines of our VLT data to measure the spherical degree of the 
largest amplitude pulsations. We apply the same technique to the previous Keck observation 
of G 29-38 (VKOO, COO) to prove that our method produces the same results as the raw 
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Figure 3. The chromatic ampHtudes of the six largest pulsation modes present in the VLT data. These chromatic amplitudes 
are created from Gaussian smoothed [a = 2.5A) spectra. The smoothing improves the signal to noise in each wavelength bin. 

chromatic amplitudes and to be able to directly compare the measurements of the two data 
sets. 



4.1 How to measure spherical degree 



We follow the procedure outlined by 



Thompson et al. 



(120041 ) for fitting each spectral line 



with the combination of a Gaussian and a Lorentzian superimposed on a linear continuum. 
We fit each spectrum only allowing the area of the Gaussian and Lorentzian to vary along 
with the slope and intercept of the linear continuum. Prior to fitting we establish the initial 
conditions of all the line parameters by fitting the average spectrum; this established the 
fuU-width-half-max values of both functions and the central wavelength of the line. 

As previously discussed, one method to present the line shape variations is to create chro- 
matic amplitudes by measuring the amplitude of each pulsation mode for each wavelength 
bin. In this case we use the fitted spectra to create the chromatic amplitudes. Figure 4 shows 
the chromatic amplitudes for a temperature and gravity model consistent with G 29-38. The 
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fits to the model spectral lines are over a limited wavelength range and do not measure 
the amplitudes between the spectral hues. As expected each spherical degree has a different 
chromatic amplitude shape. 

Alternatively, since we fit the spectra, we can reduce the problem of distinguishing I 
to quantities measured from the varying fitted parameters. Specifically, we measure the 
fractional amplitudes of the fitted EWs for the Gaussian {Gew) and the Lorentzian (Lew) 
functions at the known periods of pulsation. We also measure the fractional light amplitude 
for each pulsation mode from the fiux measured by the linear continuum of the fit at the 
central wavelength of the line (A^). These fiux amplitudes agree within the errors with those 
measured directly from the spectra in Table 1. We normalize the EW amplitudes by this 
light amplitude for each mode and plot Gew / ^\ against Ijew / ^\ for each spectral fine in 
Figure 4. 

For the data and models considered here we fit across 200 A for the H/? and H7 lines 
and 100 A for the H(5 line. All Doppler shifts of the spectral line have been removed prior 
to fitting the spectra to reduce the number of parameters in the fit, leaving the four pre- 
viously mentioned fitted parameters. Generally, the Lorentzian fits the wide wings of the 
spectral lines while the Gaussian is limited to the central core. The Gaussian represents only 
about 3 per centof the line area, however it also shows larger fractional variations than the 
Lorentzian. 



Using the limb-darkened models of DA white dwarf stars (jWeidemann fc Koesterlll984l ). 
we pulsate the models with a spherical harmonic according to ^=1,2,3 and 4. Each model 
is Gaussian smoothed by cr = 2.9 A to account for the seeing that reduces the resolution 
of the spectra, though this has only a small effect on the results. We force the fits of the 
model spectra and the observations to have the same ratio between the full-width-half- 
max of the Gaussian and the Lorentzian functions. In this way we further control the fit, 
providing consistent results between the observations and the model. We choose a model of 
T=ll,750 K and log(g)=8.25 (cgs), appropriate given the temperature and gravity measured 
by VKOO. 

Fits to model spectra show that different values of £ will produce unique Gaussian and 
Lorentzian EW amplitudes. A comparison of the chromatic amplitudes and the EW varia- 
tions in Figure 4 show the correlation between the two methods of presenting the periodic 
line shape variations. The amplitude of the Gaussian EW corresponds to the central bump 
seen in the chromatic amplitudes, while the Lorentzian EW corresponds to the wings of the 
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Figure 4. Analysis of model time series spectra for T=11,750K and log(g)=8.0 pulsated with <?=1,2,3 and 4. Each frame 
represents the fit to a different hydrogen Balmer line. Top figures are chromatic amplitudes. Bottom are Gaussian vs. Lorentzian 
EW fractional amplitudes (Gew and Ijew) normalized by the light fractional amplitude (A;^) measured for each line. See i|4.1l 
for details. 

chromatic amplitudes. With only two parameters, the EW variation plot simplifies the com- 
parison between observed pulsations and the models. In this paper we use the normalized 
EW variation plots to quantify the Keck, VLT and model line shape variations. 



4.2 Measured EW variations 

We measure the periodic changes in the normalized EWs for the six largest modes of the 
VLT observations and perform the same analysis on the Keck time series spectroscopy of 
G 29-38 (VKOO). With a consistent analysis we can compare the data sets to the models as 
well as to each other. Figures 5 and 6 show analysis of both data sets; they are consistent but 
show some offset from the model. The presented error bars are the formal errors as measured 
from sine wave fits at the stated period to the time series of the EWs and continuum fiuxes 
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(in the next section we describe monte carlo simulations that demonstrate these errors are 
reasonable.) 

Using chromatic amplitudes measured straight from the data COO discovered five 1=1 
modes and one £=2 mode. Instead, we fit the spectra to create EW variation plots for 
the 6 largest, real modes they measured in their light curve. Our results from fitting the 
spectra agree with the previous conclusions of COO. Figure 5 shows these results along with 
the model. While a small offset between the model and the Keck observations exist, the 
general trends between different spherical degree appear to be accurate. These trends are 
independent of the model temperature and gravity and thus, we can still determine C. because 
we can compare observed modes. Specifically, the mode labelled 'd' is the mode that COO 
measured as £=2 (776 s). It has larger variations in the Lorentzian and Gaussian EW than 
the other observed modes, exactly as seen in the model. 

Our results from fits to the VLT spectral series taken of G 29-38 are shown in black in 
Figure 6. Only the six modes with the largest light amplitudes are shown. The majority as- 
suredly overlap with the cluster of modes observed with the Keck data (circles). Specifically, 
our Fl and F2 are the same period as 'a' and 'b' in Figure 5, and they agre e in value within 



the e rror bars. The majority of modes in G 29-38 are expected to be £=1 ( jKleinman et al 



19981 ). and our observed cluster in the Keck and VLT observations agree with this conclusion. 
We note two interesting modes in the VLT data. One is the 681 s mode, labelled 'F6'. 
It consistently shows larger EW amplitudes in all the spectral lines, just like the previously 
known C.=2 mode ('d' in Figure 5). The other feature is the 353 s mode, labelled 'F4'. It 
consistently shows smaller EW variations, similar to what might be expected of an £=4 
mode. We note that this mode lies very near to the combination of the two primary modes 
(F1+F2) and the small combination could be influencing our results. However, simulations 
that include nearby combinations modes have no effect on the measured EW variations. 



4.3 Noise simulations 

To better understand how noise affects our measurement of spherical degree, we perform 
Monte Carlo simulations of model spectral series. Gaussian noise was introduced to both 
the overall brightness of the spectra and across the individual wavelengths of the model 
spectrum to match the noise observed in the VLT data. The error produced in individual fits 
is equivalent to the error observed in the same amplitude mode in the VLT data, confirming 
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Figure 5. Normalized EW amplitudes of each spectral line (H/3, H7 and H(5) of the Keck data as compared to the model. 
Similar to Figure |4] The model is labelled according to the spherical degree it represents (1, 2 and 4). The Keck modes are: 
a=614 s, b=818 s, c=653 s, d=776 s, e=283 s, f=430 s. The mode labelled 'd' is the one claimed by COO to have a spherical 
degree of two. 
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periods associated with each labelled mode. 



that we are closely emulating the observations. We produce and fit the H7 line of 500 
simulated spectral series (with 604 spectra in each series) and plot regions containing the 
central 68 per cent of the points for each L We show this scatter in Figure [6] for amplitudes 
of 27 mma, 10 mma and 8 mma, consistent with modes Fl, F4 and F6 respectively. 

From the simulations for the VLT noise level, we expect to be able to distinguish i* = 1, 
2 and 4 for the largest amplitude modes even if we only measured the gamma line. Given 
the simulated scatter, we would not be able to distinguish F6 from a single spectral line. 
However, the fact that we see the same trend in all three lines, increases the probability that 
we correctly identified all our modes. The simulations show that the normalized Gaussian 
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Figure 7. Results of a 500 iteration monte carlo simulation of fits to H7 of model spectra emulating the noise of the the VLT 
observations. For each t, 68 per cent of the 500 iterations lie in the area outUned. The noise free model is labelled by the number 
of its I value. Fl, F4 and F6 of the VLT measurements are included for comparison. The extent of the scatter is similar to the 
error bars of the equivalent amplitude mode, however the Gaussian and Lorentzian measurements are correlated. 

and Lorentzian EW amplitudes are correlated, but the size of the error bars are basically 
reflected by the la error bars presented in Figures 5 and 6. 



4.4 Discussion 

Due to the limb-darkening variation across the spectral lines, periodic changes in the line 
shape reveal the distribution of the pulsation across the surface of the star. The observations 
of the line shape do not explicitly agree with the model, as seen in previous chromatic 
amplitudes (COO). However, the trends between different spherical degrees do agree and as a 
result, we have measured ^ on G 29-38. The Fl, F2, F3 and F5 are all £ = 1. F6 (681 s), with 
its larger EW amplitudes, is consistent with the £=2 model and the previously identified £=2 
mode at 776 s (COO). F4 shows decreased EW amplitudes in each spectral line, indicative of 
£=4. 

We note two reasons to be cautious of F6's identification as C.=2. First, F6 does not 
show the large R^, value seen by VKOO for their £=2. Considering the error on the flux and 
velocity amplitudes, this could simply be the result of noise. Second, the 2a error bar of 
F6, as determined from the monte carlo simulations, includes the 1=1 model. The observed 
increase in EW amplitudes of one spectral line could simply be due to the scatter. However, 
since we see the same trend in all three spectral lines, it is likely that this is a real effect 
and we have correctly identified F6. 

F4 (353 s) consistently has smaller EW variations in each spectral line than every other 
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mode, similar to the models of £=4. Only one other £=4 mode has been measured on a 
DAV. In that case th e spectroscopic evidenc e was bolstered by UV observations and com- 
bination amplitudes (jThompson et al.ll2004j ). Like our observations, this £=4 on G 185-32 
shows reduced G aussian and Lorentzian EW amplitudes as compared to the other modes 
(JThompsonI l2006l ) . The similarity between our observations of F4 and this established £=4 
mode gives strength to the claim that F4 is £=4. 

Another indicatio n of higher d, i s the p resence of large am plitude har monics as compared 



to the parent modes. 



Yeates et al. 



Wu 



(J200ll ) along with the am- 



(120051 ) use the theories of 
plitudes of combination modes to show that G 185-32 does indeed have a higher I mode, 
though their results better agree with £=3. Here, F4 has harmonic modes despite its low 
amplitude, while the larger amplitude mode, F3, shows no obvious harmonics. More specif- 
ically. Re for 2F4 is significantly larger than 2F1 (the harmonic of an £ =1 mode) . However, 
2F4's Re value is closer to what we would expect for an i=2 harmonic (jWull200ll ). Analysis 
of the combination modes from the VLT data set is somewhat suspect since they all have 
low amplitudes and are easily influenced by noise. Plus, relative amplitudes of combination 
modes are influenced by azimuthal order. W ith such a short data set we cannot measure the 
azimuthal splitting (jKleinman et al.l 119981 ) nor its effects on the combination amplitudes. 
While the combinations enforce the idea of higher £ for F4, they leave some question to its 
actual identity. 

Part of the ambigu ity in the identificatio n of F4 may be caused by our entirely linear 
treatment of pulsation. Ilsing &: Koesterl (l200ll ) show that pulsation modes can be distributed 
across the star in a manner significantly different from a simple spherical harmonic. From 
their simulations, they find that an £=3 with large amplitudes can have chromatic amplitudes 
more like i=l or 4 than £=3. Higher C. modes must have larger amplitudes in order to be 
observed, and are more likely to suffer from these nonlinear effects. This could be the reason 
we observe line shape variations similar to £=4 modes, but not £=3; nonlinear mixing distorts 
the mode identification. In addition to our F4, these nonline arities may b e respo nsible for 



Kotak et al 



(|200J). Luckily, 



the unusual chromatic amplitudes seen on G 117-B15A by 

knowing that the mode is high C. is enough to improve asteroseismological modeling. Given 

these considerations, we cautiously denote F4 as d =4 or 3. 

Because F4 lies near to the combination of Fl and F2, it is wo rth considering the 



possibility that 



Wu fc Goldreich 



'4 is a naturally damped mode excited by resonance (JDziembowski 



1982; 



20011 ). However the three-mode coupling angular selection rules show that 
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Table 2. Mo de identifications for G 29-38 using time series spectroscopy. VLT refers to those modes identified in this paper. 
COO referes to lClemens et ahl l|2000h and K02 referes to Kotak etftl. tSQOi). See ^Xil for more details on the 353s mode. 



period 


source 


£ 


937 


VLT 


1 


920 


K02 


2 


835 


VLT 


1 


815 


COO.VLT 


1 


776 


COO 


2 


681 


VLT 


2 


655 


COO 


1 


614 


COO.VLT 


1 


431 


COO 


1 


353 


VLT 


4 or 3 


284 


COO 


1 



because Fl and F2 are both £=1, F4 could only be excited by this mechanism if it were £=2 
( iWu &: Goldreiclj|200ll ). Given our measurements of the line shape variations, resonance is 
unlikely in this instance. 

In Table 2 we give a summary of all mode identifications from time series spectroscopy 
of G 29-38. We have eleve n identificati o ns, si x from the Keck data made by COO, one from 

(120021 ). and four new identification (plus two confir- 



the Keck data inferred by 



Kotak et al. 



mations) from our VLT data. Note that we have identified modes to be the same if their 
fre quencies matched to w ithin ~ 10 /xHz, which is roughly the rotational splitting inferred 



by 



Kleinman et al. 



mm . 



5 CONCLUSIONS 

Time series spectroscopy has proven to be a useful tool to measure white dwarf pulsation. 
For G 29-38, it has revealed the physical motion of the pulsation through Doppler shifts and 
a through the variations of line shape. The three sets of spectroscopy on G 29-38, the VLT 
time series present ed here, the Keck time series (VKOO and COO), and the high- resolution 



Keck observations (jThompson et al 



20031 ). have proven that this technique gives consistent 



results and allows us to measure aspects of the star's pulsation that cannot be obtained 
through photometry. 



The convective driving theory (JGoldreich fc Wu 



1999 



Brickhill 



19831 ) is now the favored 



theory to produce pulsat i ons on DAVs. The velocity modes that ar e measured here, i n VKO O 



and in 



Thompson et al.l (120031) show amplitudes as predicted by IWu fc GoldreichI (1l999l ) 



However, we still lack the appropriate data to test the predictions at higher frequencies. 
R^ should increase at shorter periods, but G 29-38 sadly does not tend to drive pulsation 
modes shorter than 400 s as strongly as it does 800 s modes, making it hard to measure 
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the trend without extensive time series spectroscopy. Further spectro scopic observations are 
also warranted because of the velocity combinations (our FH-F2 and [Thompson et al.ll2003h 
not anticipated by any theory. Further observations will determine the general properties of 
these velocity combinations whether they represent nonlinearities in the velocity curve or 
are independent pulsations possibly driven by resonance. 

The major achievement of time series spectroscopy is the measurement of £, necessary 
for asteroseismology. This technique promises to reveal internal properties s uch as the hy- 



Bradley 



1998; 



drogen layer mass, helium la yer mass, temperature, mass etc. (see for example 
Castanheira fc Kepleiil2008l ). Modeling white dwarf interiors requires knowledge of £, oth- 
erwise more than one predicted mode could easily fit each observed period, and more than 
one solution exists for the star. Without any other information, £=1 is generally assumed 
for all modes, except for those not easily fit by £=1. As we have shown for G 29-38, that 
assumption is frequently incorrect; 4 of the 11 measured modes are d^l. Past attempts to 
model G 29-3 8 resulted in an abundance of possible solutions that could fit the complex set 
of pulsations (jKleinmanl Il995l ) . Now that we have constrained £ for G 29-28, asteroseismic 
modeling is more li kely to succeed. Thou gh beyond the scope of this paper, by combining 
known photometry (JKleinman et al.lll998l ) with the d, values presented here, we expect that 
asteroseismological modeling will finally be able to determine the internal structure of the 
notorious G 29-38. 
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